Two-proton correlations from 158 AGeV Pb+Pb central collisions 
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The two-proton correlation function at midrapidity from Pb+Pb central collisions at 158 
AGeV has been measured by the NA49 experiment. The results are compared to model pre- 
dictions from static thermal Gaussian proton source distributions and transport models rqmd 
and VENUS. An effective proton source size is determined by minimizing \ /ndf between the 
correlation functions of the data and those calculated for the Gaussian sources, yielding <r e g = 
3.85 ± 0.15(stat.)^g'2g(syst.) fm. Both the rqmd and the venus model are consistent with the data 
within the error in the correlation peak region. 
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Nuclear matter at high energy density has been exten- 
sively studied through high energy heavy ion collisions 
(for recent developments, see fl). It is hoped that these 
collisions will create a deconfined state, the quark-gluon 
plasma. These studies have shown that nucleus-nucleus 
collisions are not mere superpositions of nucleon-nucleon 
collisions. In Pb+Pb central collisions at the CERN SPS, 
the nucleons are transported further into the midrapidity 
region than in S+S central collisions and p+p interac- 
tions [^] ^) , consistent with the picture that the incoming 
nucleons in nucleus-nucleus collisions have undergone a 
significant number of scatterings. This results in a high 
rapidity density of net-baryons at midrapidity The 
spatial baryon density plays an important role in the dy- 
namical evolution of heavy ion collisions |(| . An essential 
ingredient for establishing the spatial baryon density is 
the space-time extent of the baryon source at freezeout, 
which can be inferred from two-proton correlation func- 
tions §,| . 

Two-proton correlations at low relative momentum 
are due to final state interactions (the attractive strong 
and the repulsive Coulomb interactions) and Fermi-Dirac 
quantum statistics . The correlation function is zero 
at <7i, w = = and peaks at g; nv sw 20 MeV/c, 

where is the difference of the proton 4-momenta, and 
Qinv is the momentum magnitude of one proton in the rest 
frame of the pair. The peak height is inversely related to 
the space-time extent of the proton source . 

The space-time extent of the pion source has been stud- 
ied extensively via two-pion interferometry in heavy ion 
collisions [p|~pd|. However, measurements of the space- 
time extent of the proton source are rare, especially for 
high energy nucleus-nucleus collisions fll2fl . In this letter, 
we report the first measurement of the two-proton cor- 
relation function in the midrapidity region from Pb+Pb 
central collisions at 158 AGeV, by the NA49 collabora- 
tion at the CERN SPS. 

NA49 is a large acceptance hadron spectrometer [ ^3| . 
The main detectors are four large time-projection cham- 
bers (TPCs). Two TPCs (VTPCs) are placed along the 
beam axis inside two dipole magnets, which have a max- 
imum integrated field strength of 9 Tesla-meters. The 
other TPCs (MTPCs) are placed downstream of the mag- 
nets on either side of the beam axis. Behind the MTPCs 
are two time-of-flight (TOF) walls, covering smaller ac- 
ceptance than the TPCs. A beam of 208 Pb struck a Pb 
target of thickness 224 mg/cm 2 , placed in front of the 
first VTPC. A zero-degree calorimeter, located further 
downstream on the axis of the deflected beam, measures 
the kinetic energy of the projectile spectators. By requir- 
ing less than 6 TeV energy measured in the calorimeter, 
the 5% most central events were selected. These events 
correspond approximately to collisions with impact pa- 
rameter 6 < 3.3 fm. 

Two independent analyses of the two-proton correla- 
tion function were carried out: 



1. dE/dx analysis: Track segments reconstructed in 
the MTPCs were matched to track segments from 
at least one of the VTPCs. Particle momenta were 
determined by track curvature in the VTPCs. The 
particle identification was performed by measur- 
ing the specific ionization (dE/dx) deposited by a 
charged particle in the gas of the MTPC, measured 
in the region of the "relativistic rise" . The mean 
(dE/dx) of a particle was estimated using the trun- 
cated mean technique |l4j . A relative (dE / dx) res- 
olution of it ~ 5% was achieved in this analysis. 
Since the proton (dE/dx) is approximately lcr and 
3d below those of the kaons and the pions with the 
same momentum, protons cannot be uniquely iden- 
tified. Instead, particles with at least 70% proba- 
bility to be a proton were included in the analysis. 

2. TOF analysis: Only tracks reconstructed in the 
MTPCs were used. The momentum of a particle 
was determined by an iterative procedure in which 
the trajectory fitted to the measured points was 
projected upstream through the magnetic field, as- 
suming that it originated at the primary interaction 
vertex. The particle identification was performed 
by combining the momentum measurement with 
the velocity information from the TOF walls and 
the (dE/dx) from the MTPCs. The typical TOF 
resolution was measured to be 60 ps. 
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FIG. 1. NA49 proton acceptance in transverse momentum 
Pt versus rapidity y. The acceptance values for protons used 
in the dE/dx analysis are shown in the box plot between the 
dashed lines, and for the protons used in the TOF analysis 
are shown in the contours. 

Fig. [I] shows the proton acceptances from the two anal- 
yses, which are in different hemispheres with respect 
to midrapidity (y c . m . = 2.9). The two-proton correla- 
tion function was analysed using protons from the ra- 
pidity range 2.9 < y < 3.4 in the dE/dx analysis, and 
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2.4 < y < 2.9 in the TOF analysis. Both analyses used 
protons up to a transverse momentum px = 2 GeV/c. 

The two-proton correlation function was calculated as 
the ratio of the <7i nv distribution of true proton pairs to 
that of mixed-event pairs with protons taken from dif- 
ferent events. The number of mixed-event pairs is large 
enough that the statistical error on the correlation func- 
tion is dominated by the statistical uncertainty in the 
number of true pairs. To eliminate the effect of close 
pair reconstruction inefficiency, a cut of 2 cm was ap- 
plied on the pair distance at the mid-plane of the MTPC 
for both true and mixed-event pairs |0| . 

The proton sample is contaminated by protons from 
weak decays (A including feeddowns, and E + ) which are 
incorrectly reconstructed as primary vertex tracks. From 
the measured single particle distributions dJOl as weu 
as model calculations of RQMD and VENUS, we estimate 
the contamination to be 25 ± 10% (with weak pt de- 
pendence). This results in 44^ 6 % of the proton pairs 
having at least one proton from weak decays, which is 
assumed to have no correlation with the protons pro- 
duced in the primary interaction [ jl7| . Hence, protons 
from weak decays only reduce the correlation strength, 
without changing the shape of the two-proton correlation 
function at low qi nv . Additional contamination is present 
in the dE/dx analysis from kaons on the lower tail of their 
(dE/dx) distribution merging with the region where par- 
ticles have at least 70% probability to be a proton. This 
results in 25% K + p pairs and fewer than 2% K + K + pairs 
in the proton pair sample, independent of ow- K + p 
pairs with < k < 100 MeV/c (k is the momentum of 
the K + or proton in the rest frame of the pair), which 
spans the region affected by Coulomb repulsion for the 
typical source created at midrapidity in Pb+Pb central 
collisions, affect the two-proton correlation function in 
the range 160 < q lmr < 460 MeV/c. Only K + p pairs 
with 180 < k < 240 MeV/c contribute to the two-proton 
correlation function in the range < (jw < 50 MeV/c. 
Therefore, the K + p contamination in the proton pair 
sample shows no structure in qi nv below 160 MeV/c, but 
reduces the correlation strength by 25%. The contribu- 
tion from the K + K + pairs in the two-proton correlation 
function is negligible. 

The correlation functions obtained from the two anal- 
yses can be directly compared because of the nearly sym- 
metric acceptances used. The corrected correlation func- 
tions are consistent. In the results reported below, the 
corrected true pairs and the respective mixed-event pairs 
from the two analyses are combined. The combined sam- 
ple has about 10 5 pairs with q inv < 120 MeV/c, 75% of 
which are from the TOF analysis. The gw distribu- 
tions of the true and the mixed-event pairs are shown 
in the top panel of Fig. |2| The number of mixed-event 
pairs is normalised to that of true pairs in the range 
qinv > 500 MeV/c. The resulting correlation function 
C raw is shown in the lower panel of Fig. ||, and tabulated 



in Table ffl. The peak at q luv ~ 20 MeV/c has amplitude 
1.14 ±0.04. A statistically significant structure is seen in 
the correlation function at qi nv « 70 MeV/c. The struc- 
ture is present in both the dE/dx and TOF data sets. 
See below for further discussion of the structure. 
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FIG. 2. Upper panel: (jw distributions of true proton pairs 
(points) and mixed-event proton pairs (histogram). Lower 
panel: the measured two-proton correlation function. In both 
panels, the points are plotted at the IllCclIl Qinv of each bin. 
The bin size is 6 MeV/c. The errors shown are statistical 
only. The contamination from weak decay protons and the 
finite qi nv resolution are not corrected. 

The correlation function was further corrected for the 
contamination from weak decay protons and the finite 
qinv resolution. The former is corrected via C = (C law — 
0.44)/(l - 0.44). The momentum resolution is 0.25% 
at a momentum about 10 GeV/c; the yi nv resolution is 
2.5 MeV/c in the range q mv < 120 MeV/c |||. The ef- 
fect of the (?i nv resolution was studied using the RQMD 
and VENUS models, and was found to be significant only 
for the first two data points of the measured correlation 
function. The final corrected correlation function is plot- 
ted in Fig. U as filled points, and tabulated in Table || 

In order to assess the proton freeze-out conditions, we 
compare the measured two-proton correlation function to 
theoretical calculations. Given the proton phase space 
density distribution, the two-proton correlation function 
can be calculated by the Koonin-Pratt Formalism 0,^8) . 
The formalism uses the phase-shift method, which incor- 
porates the Coulomb and the strong interactions between 



3 



the protons. We use two types of proton freeze-out dis- 
tributions: 

(I) Gaussian sources of widths a x<y ,a z and at for the 
space and time coordinates of protons in the source rest 
frame, and a Boltzmann distribution with temperature 
T for the proton momentum distribution. No correla- 
tion between space-time and momentum is present. The 
following combinations of parameters are used: a x y — 
a z =a\ a t = and a; T — 120 MeV (as derived in @), 
300 MeV (measured inverse slope of proton transverse 
mass spectrum and 70 MeV (inverse slope observed 
at low energy, as an extreme). 

(II) Protons generated for Pb+Pb central collisions 
(b < 3.3 fm) at 158 AGeV by two microscopic trans- 
port models: the rqmd model (version 2.3) jL9| and the 
VENUS model (version 4.12) pQ[. Both models describe 
a variety of experimental data on single particle distri- 
butions reasonably well. Protons at freeze-out have cor- 
relations between space-time and momentum intrinsic to 
the dynamical evolution in the models. Neither model 
includes correlations due to quantum statistics and final 
state interactions. Particle weak decays are not included 
in the models. 
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FIG. 3. The two-proton correlation function after correc- 
tions for the 44% contamination due to weak decay protons 
and the finite qinv resolution (points), compared to calcu- 
lations for a Gaussian source (dotted histogram), and for 
freeze-out protons from rqmd v2.3 (solid histogram) and 
VENUS v4.12 (dashed histogram). The errors shown on the 
data points are statistical only. The dot-dashed line indicates 
the Qinv value (48 MeV/c) below which the comparisons be- 
tween the data and the calculations are performed. 

Only protons in the experimental acceptance are used 
to calculate the two-proton correlation functions, the re- 
sults of which are shown in Fig. ||. The correlation func- 
tion for the Gaussian source with a x>y = a z = 3.8 fm, 
at = fm and T = 120 MeV, which describes well the 



experimental data, is also plotted. 

We use % 2 /ndf, the normalised mean square of the 
point-to-point difference between the data and the cal- 
culation in the range < 48 MeV/c (i.e., 8 data 
points or ndf = 8), to quantify how well the calcu- 
lations agree with the data. We characterize the ef- 
fective size of the model-generated proton source by 
CTcff = $/aAx ■ a Ay ■ aAz/V2, where aAx,aA y and aAz 
are the Gaussian widths fitted to the distributions in 
Ax, Ay and Az, the distance between the protons of close 
pairs with q\ nv < 48 MeV/c. The distance is evaluated in 
the pair rest frame (since the correlation function is stud- 
ied as a function of q- mv ) at the time when the later par- 
ticle freezes out [pi). Repectively for RQMD and VENUS, 
the x 2 / n df values are 1.53 and 1.38, corresponding to 
a probability of 20% and 14% that the model distribu- 
tions are consistent with the data; the fitted Gaussian 
widths are (a&. x , a a v , aAz) — (5.91,6.00,6.83) fm and 
(4.57,4.57,6.08) fm, where z is the longitudinal coordi- 
nate; consequently the effective sizes are a c g = 4.41 fm 
and 3.55 fm. 



3 





1 ' 


' l 


i | i i i i | i 








•— o 


T= 120 MeV 








A A 


T = 300 MeV 








* in- 


T= 70 MeV 








filled 


o x =o y =o z =a, 








open 


a x =O y=°z. CI t= 






o\ 






O 








\a' o'' 


























b 


y i 










5 vol 


m 


RQMD V2.3 








□ 


Venus v4.12 








o-c T=120MeV,o x = 


=c y =4a z ,o,=0 




i , 




o o T=120MeV,o x = 


=o y =aJ4,c,=0 

'I''' 




3 


4 


5 6 


7 i 



a „ (fm) in pair rest frame 



FIG. 4. The x 2 / n df values as function of the effective 
source size a e g for various model calculations with respect 
to the measured correlation function. The model calculations 
are for Gaussian sources (circles, triangles, and stars, both 
filled-in and open) and for freeze-out protons generated by 
rqmd v2.3 (filled square) and VENUS v4.12 (open square). 
The curves connect the various sets of points to guide the 
eye. The errors on cr e ff are negligible for all the points. 

In Fig. [|, we study the \ 2 1 n df as a function of fj c ff • The 
% 2 /ndf values for all three models follow roughly the same 
solid curve, drawn through the points for the Gaussian 
sources with T = 120 MeV to guide the eye ]22|| . From 
the minimum % 2 /ndf (=0.86) point and the points where 
% 2 /ndf has increased by 0.125 (note ndf = 8), we extract 
Coff — (3.85 ± 0.15) fm, where 0.15 fm is the statistical 
error p3fl . By applying a correction to the measured 
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correlation function using a proton pair contamination of 
28% and 58%, we obtain a systematic error of ±0.15 fm 
on <7 c ff. We note that <7 e ff = 3.85 fm corresponds to a 
uniform density hard sphere of radius ^/ba c g = 8.6 fm, 
which is larger than the size of the colliding Pb nuclei. 

We have also studied Gaussian sources with extreme 
shapes: oblate a x>y — Aa z and prolate a x ^ y — a z /4 (both 
with at = 0). The corresponding x 2 /ndf versus a c g 
curves are shown in Fig. [I| as dashed and dotted lines, 
respectively. They do not fall along the curve for the 
isotropic Gaussian sources, implying that the two-proton 
correlation function has certain sensitivity to the shape of 
the proton source. The data do not favour the extreme 
Gaussian shapes. In principle, multi-dimensional two- 
proton correlation functions could reveal the shape of the 
proton source |Q,|l8|. However, this requires greater sta- 
tistical precision than is available with the present data 
set. 

By minimizing the x 2 /ndf between the correlation 
functions for the various Gaussian sources with respect to 
those for RQMD and VENUS, we can extract the effective 
proton source size for the models. The results are consis- 
tent with the calculated effective sizes from the models 
themselves. 

We note that there is no simple relation between a e g 
and the proton source in the models. From single proton 
distributions at freeze-out, we obtain the following Gaus- 
sian widths in the source rest frame: (a x , a y , a Zl at) — 
(7.6,7.7,6.4,7.0) fm for rqmd and (3.6,3.6,4.3,1.9) fm 
for VENUS, respectively. The two-proton correlation func- 
tion, therefore, appears to measure a smaller region of the 
source due to space-time-momentum correlation. The ef- 
fect is more dramatic in RQMD than in VENUS, which is 
consistent with the expectation that more secondary par- 
ticle interactions in rqmd result in a stronger correlation 
between space-time and momentum of freeze-out pro- 
tons. The fact that the x 2 / n df values versus a e g for the 
models lie on the curve obtained for the Gaussian sources, 
in which no space-time-momentum correlation is present, 
suggests that the effect of the space-time- momentum cor- 
relation is small in a e g. 

Now we come back to the structure observed at qi nv w 
70 MeV/c in the correlation function. Many systemat- 
ics have been checked, including (1) pt dependence, (2) 
stricter acceptance cut, (3) larger cut on two track dis- 
tance, (4) varying cuts to allow different kaon contamina- 
tions in the TOF data set p7| , (5) allowing only one par- 
ticle from the same event in event mixing, (6) residual ef- 
fect in single particle spectra due to two-particle correla- 
tion, (7) misreconstructed momenta for A decay protons, 
and (8) reflection of Ap strong interactions in two-proton 
correlation with one proton from A decays |24]. None 
of the checked systematics can account for the structure. 
Nevertheless, if we assume that the structure has an os- 
cillatory behavior, cos(aq lnv ) / q- lnv (where a is a constant) , 
and also contributes to the correlation peak region, then 



we estimate a systematic error of t-oxo ml m the ex- 
tracted <7 e ff, in addition to the ±0.15 fm estimated above. 

It is, however, possible that the structure is due to 
underlying physics. For instance, such structure may re- 
sult from a breakdown of the basic assumptions under- 
lying the Koonin-Pratt Formalism ])25||: the smoothness 



assumption |26|, the independent emission assumption, 
and/or the q- lnv independence of the two-proton source. 
A hard edge in the two-proton source distribution may 
also result in similar structures, however, it does not re- 
produce the observed sign and amplitude, unless the in- 
dependent emission assumption is violated. See also the 
discussion in Ref. Q . 

Finally, we comment on our two-proton correlation 
function in the context of other measurements. The pion 
source size measured by interferometry increases with the 
pion multiplicity |p"i"|| , which increases steadily with bom- 
barding energy in similar colliding systems f27|| . Due to 
the large pion-nucleon cross-section, one would expect 
that protons and pions freeze-out under similar condi- 
tions, therefore, the proton source size would increase 
with bombarding energy as well. However, our measure- 
ment, in conjunction with preliminary results obtained 
at GSI |^8| and AGS p^] energies, shows that the peak 
height is rather insensitive to the bombarding energy. 
This implies that the effective sizes of the freeze-out pro- 
ton sources are similar in heavy ion collisions over a wide 
energy range. 

In summary, NA49 has measured the two-proton cor- 
relation function at midrapidity from Pb±Pb central col- 
lisions at 158 AGeV. The peak height of the two-proton 
correlation function is 1.24±0.07 after corrections. From 
comparison between the data and the calculations for 
static thermal Gaussian sources, we extract an effective 
proton source size <7 e ff = 3.85±0.15(stat.)^Q2g(syst.) fm. 
Within the error, the RQMD model (a c g = 4.41 fm) and 
the VENUS model (a c g = 3.55 fm) are consistent with 
the data. Due to the space-time-momentum correlation, 
the two-proton correlation function is sensitive only to a 
limited region of the proton source. Our measurement to- 
gether with the measurements at lower energies suggest 
a very weak dependence of the two-proton correlation 
function on bombarding energy. The observed structure 
at (fr nv 70 MeV/c is not understood. 
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TABLE I. The two-proton correlation function. (?i n v): 
mean Qjnv 

in MeV/c; C ra wi uncorrected correlation function; 
C: correlation function corrected for 44% contamination in 
the proton pair sample due to weak decay protons and for the 
finite q, m resolution. All quoted errors are statistical only. 



(?inv) 


Craw 


C 


4.8 


0.500 ± 0.224 


0.084 ±0.314 


9.7 


0.949 ± 0.092 


0.892 ±0.161 


15.4 


1.049 ±0.055 


1.087 ± 0.099 


21.3 


1.136 ±0.041 


1.244 ± 0.074 


27.2 


1.051 ±0.031 


1.090 ± 0.056 


33.2 


1.033 ±0.026 


1.058 ± 0.046 


39.1 


0.965 ±0.021 


0.937 ± 0.038 


45.1 


0.956 ±0.019 


0.921 ± 0.033 


51.1 


0.966 ±0.017 


0.940 ± 0.030 


57.1 


0.987 ±0.015 


0.977 ± 0.028 


63.1 


1.007 ±0.014 


1.012 ±0.026 


69.1 


1.027 ±0.013 


1.048 ± 0.024 


75.1 


0.998 ±0.012 


0.996 ± 0.022 


81.1 


0.976 ±0.012 


0.958 ± 0.021 


87.1 


0.974 ±0.011 


0.954 ±0.019 


93.0 


0.993 ±0.010 


0.988 ±0.019 


99.0 


1.005 ±0.010 


1.010 ±0.018 


105.0 


0.998 ±0.010 


0.996 ±0.017 


111.0 


0.990 ±0.009 


0.982 ±0.017 


117.0 


1.001 ±0.009 


1.002 ±0.016 
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